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Abstract 

We discuss the angular distribution of dilepton pairs in p+A collisions at LHC 
at low momentum transfer and transverse momentum. This includes helicity 
amplitudes and angular coefficients as well as a look at the Lam-Tung sum 
rule. We use leading twist NLO calculations and higher twist corrections by 
double scattering in large nuclei. 


1. INTRODUCTION 

The Drell Yan cross section can be described by a contraction of the lepton and hadron tensors 
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Here we parametrize the lepton pair by the invariant mass Q of the virtual photon and its transverse 
momentum q± and rapidity y in the center of mass frame of the colliding hadrons. In addition we give 
the direction of one of the leptons, say the positively charged one, in a photon rest frame using polar and 
azimuthal angles <i> and 9: dQ = dxl> d cos 9. Now the cross section can be understood in terms of four 
helicity amplitudes [|l], |j, [J] 
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These are defined as contractions W a y = e M (a)W fil 'e*(a') of the hadron tensor with polarization 
vectors of the virtual photon for polarizations a = 0, ±1. Only four out of all possible contractions are 
independent, the others can be related by symmetries of the hadron tensor. The usual choice is to pick 
the longitudinal Wl = lUo.o- the helicity flip Wa = (Wi ,o + Wq,\)/\/2 and the double helicity flip 
amplitude Waa = W\-\ together with the trace IUtl = Wr + Wl /2 = — Wj //2 as a basis. Note that 
integration over the angles 9 and 4> leaves only contributions from the trace 
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On the other hand, if we are only interested in relative angular distributions — i.e. in the ratio 
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— we can make use of angular coefficients. Two different sets can be found in the literature [|4]. One set 
consists of the coefficients 
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the other one is defined by 
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Fig. 1: The Collins-Soper frame: the 2 -axis cuts the angle between Pi and —P 2 into halves (the half angle is called the 
Collins-Soper angle 7cs) while the x-axis is perpendicular to Pi and P 2 . The direction of one lepton momentum li can then 
be given by the angles 8 and rj> 


The helicity amplitudes are frame dependent. In principle we allow all frames where the photon 
is at rest, i.e. q 1 ' = (Q. 0,0,0). However there are some frames with particular properties studied in [ 2 ]. 
Here we only use the Collins-Soper (CS) frame. It is characterized by two properties. First the y-axis is 
perpendicular to the plane spanned by the two hadron momenta Pi and P 2 (which are no longer collinear 
in a photon rest frame as long as q_ / 0, what is true in our kinematic domain) and second the 7-axis 
cuts the angle between Pi and —P 2 into two equal halves, see Fig. |l|. 

2, LEADING TWIST 

The hadron tensor for the Drell Yan process in a leading twist (twist-2) calculation is given by the well 
known factorization formula as a convolution of two parton distributions with a perturbative parton cross 
section. We are interested here in the kinematic region characterized by intermediate photon mass Q and 
intermediate transverse momentum q± ~ Q of a few GeV. For these values the dominant contribution 
to twist-2 is given by the next-to-leading order (NLO) perturbative diagrams like the one in Fig. (left) 
and we can safely omit logarithmic corrections of type In 2 Q 2 /q]_. It has been shown that a leading twist 
calculation up to NLO respects the so called Lam Tung sum rule Wl = 2Waa PI- In terms of angular 
coefficients this can be rewritten as A$ = .4 2 or 2u = 1 — A. Furthermore the spin flip amplitude Wa 
has to vanish for a symmetric colliding system like p + p. For p + A we expect small contributions for 
Wa due to lost isospin symmetry and nuclear corrections to the parton distributions. Results for p + p at 
\/S = 5.5 TeV have already been presented elsewhere [j6j]. 

3. NUCLEAR ENHANCED TWIST-4 

For large nuclei corrections to the leading twist calculation, induced by multiple scattering, play an 
important role. The formalism how to take into account these nuclear enhanced higher twist contributions 
was worked out by Luo, Qiu and Sterman [|7j. 8. |9j]. The leading nuclear corrections (twist-4 or double 




Fig. 2: Examples for diagrams contributing to the dilepton production in hadron (h) nucleus (A) scattering at twist-2 (left) and 
twist-4 (right) level 
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Fig. 3: Rapidity dependence of ALWtl (left) and N^Wa (right) forp + Pb at -\/S = 8.8 TeV, Q = 5 GeV and q± = 4 GeV: 
twist-2 NLO (short dashed), double hard with EKS98 modifications (long dashed), soft hard (dash dotted, scaled up by a factor 
of 10) and the sum of all contributions (bold solid line). The double hard contribution calculated without EKS98 modifications 
is also shown (dotted line). Note that the incoming nucleus has positive rapidity. 


scattering) have already been calculated for some observables. For Drell Yan this was first done by Guo 
[il(]] and later generalized to the Drell Yan angular distribution | I. Tl|]. 

Fig. (right) shows an example for a diagram contributing at twist-4 level. Now two partons a, 
b from the nucleus and one (c) from the single hadron are involved. As long as q± ~ Q twist-4 is 
dominated by two different contributions. The double hard (DH) process where each parton from the 
nucleus has a finite momentum fraction and the soft hard (SH) process where one parton has vanishing 
momentum fraction. The factorization formulas are given by 







- £ r c5 H (Ta, X b ) H^ +c (q, X a , X b , X c ) fc(x c ), 

Xq 

—— D Xa! x b {Q, x c)T a b (x a ) H^ b+C {q, x a , x b = 
x c 


0, x c ) f c {x c ) 


(V) 

( 8 ) 


for double hard and soft hard scattering respectively § 1, 0. D Xa .x b (q- I x c) is a second order dif¬ 
ferential operator in x a and x b . T^ H and T^ H are new matrix elements of twist-4 which encode non 
perturbative correlations between the partons a and b. Since we are still missing solid experimental 
information about these new quantities, they are usually modeled in a simple way through parton distri¬ 
butions. We use T^ b n (x a ,x b ) = CA 4 / 3 f a (x a )f b (x b ) and T^(x a ) = A 2 A^ 3 f a (x a ) where C and A 2 
are normalization constants. The key feature, their nuclear enhancement, is their scaling with the nuclear 
size. 

It has been shown in [|4j. [TT]] that the DH process shows a trivial angular pattern in the sense that it is 
similar to the lowest order simple annihilation of on shell quarks with the only difference that one of the 
quarks now carries finite transverse momentum q±. In this spirit it is no surprise that the DH contribution 
respects the Lam Tung relation. On the other hand SH scattering is more complicated and violates the 
Lam Tung sum rule. Also we expect that the spin flip amplitude VLy can receive large contributions from 
the twist-4 calculation. 


4. NUMERICAL RESULTS 

In this section we present some numerical results obtained for proton lead collisions at the LHC en¬ 
ergy \/S = 8.8 TeV. Results for RHIC energies can be found elsewhere [|4|. We use CTEQ5L parton 
distributions [12] combined with EKS98 [13] nuclear nuclear modifications both for the nuclear par- 
ton distributions and for the models of the twist-4 matrix elements. In some plots we also give results 
for double hard contributions where the nuclear modification was omitted. Since we do not know any- 
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Fig. 4: The angular coefficients A (long dashed), y, (short dashed) and v (dot dashed) forp + Pb at = 8.8 TeV and Q = 5 
GeV as functions of rapidity (left, at q± = 4 GeV) and transverse momentum (right, at y = 0). Thin lines represent pure 
twist-2 calculations, thick lines show the results including twist-4. The solid line gives the violation of the Lam Tung relation 
2v — (1 — A) scaled up by a factor of 10. 
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thing about the correct .x-dcpcndcncc of the higher twist matrix elements this gives an impression about 
the error we may at least assume for the higher twist calculation. The normalization constants for the 
twist-4 matrix elements are chosen to be A 2 = 0.01 GeV 2 and C = 0.005 GeV 2 . In order to enable 
convenient comparison with cross sections we show all helicity amplitudes multiplied with the prefactor 
N a = a 2 / (647T 3 SQ 2 ) from Eq. (|2j). 

In Fig. H we give results for the helicity amplitudes Wtl and Wa as functions of rapidity. We 
observe that DH scattering gives a large contribution at negative rapidities (the direction of the proton) 
which can easily balance the suppression of the twist-2 contribution by shadowing. Soft hard scattering 
is strongly suppressed at these energies. Wj, (not shown) qualitatively has the same behavior as Wtl- 
The helicity flip amplitude Wa picks up only a small contribution from twist-2 and is entirely dominated 
by double hard scattering as already expected. This would be a good observable to pin down nuclear 
effects in p + A collision. Fig. [t] gives the full set of angular coefficients A, // and v as functions of y 
and q±. Both the twist-2 results and the modifications by twist-4 are given. The violation of the Fam 
Tung sum rule is numerically almost negligible since the soft hard contribution is so small. Note here 
that earlier experiments have discovered a large violation of the Fam Tung relation in 7r + A collisions 
[0]- This issue is still not fully resolved. 

p + A collisions offer the unique opportunity to study nuclear effects directly via the rapidity 
dependence of observables. Helicity amplitudes and angular coefficients can help to pin down the role 
of nuclear enhanced higher twist. Particularly promising are the helicity flip amplitude Wa and the 
coefficient // which both vanish for p + p. The usage of several different species of nuclei is advisable in 
order to address the important question of A scaling. 
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